A numerical examination of an evolving black string horizon 
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We use the numerical solution describing the evolution of a perturbed black string presented in 
to elucidate the intrinsic behavior of the horizon. It is found that by the end of the simulation, the 
affine parameter on the horizon has become very large and the expansion and shear of the horizon 
in turn very small. This suggests the possibility that the horizon might pinch off in infinite affine 
parameter. 



I. INTRODUCTION 



While in four-dimensional spacetimes black holes are 
stable, it has been shown by Gregory and Lafiamme 2] 
that at least some of their higher dimensional analogs, 
in particular black strings, are linearly unstable. As op- 
posed to the situation in four-dimensions, at the linear 
level perturbations of these black strings grow exponen- 
tially above some critical length of the string. Although 
this picture has been well understood through a variety 
of analysis 0, H S S S 0, 

it is not yet known what 
an unstable black string evolves to. In Q, using the lin- 
ear instability of the string coupled to entropy consider- 
ations, it was argued that an unstable black string would 
"pinch off" and form a set of black holes which would 
later merge. This scenario, at the classical level, must 
give rise to a violation of cosmic censorship. However, 
this possibility was made doubtful by more recent work 
of Horowitz and Maeda (HM) There, a theorem is 
proved using the properties of expansion, shear and affine 
parameter of the geodesic generators of the black string 
horizon. This theorem asserts that a pinch off in finite 
affine parameter could only occur if a singularity formed 
outside the horizon. Under the assumption that this, 
perhaps too drastic scenario, is unlikely, the only other 
option left open for the event horizon pinching-off is that 
it does so in infinite affine parameter. This possibility, 
although mathematically possible, was argued against by 
HM. Their argument is based on one main observation 
together with a plausibility argument. The observation 
is that although their theorem does indeed allow for a 
pinch-off in infinite time, the rate, with respect to affine 
parameter, at which circles at r — const can decrease 
in size is extremely slow. HM then argue that if the 
system decays so slowly, its dynamics must go through 
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a series of essentially static non- uniform strings, and it 
then might just stay at one of these solutions. Therefore, 
instead of the pinch-off scenario, HM conjectured that 
the black string would evolve into a black string that is 
static but not uniform in the z direction: a static wig- 
gly black string. Families of static wiggly black strings 
were found numerically by Wiseman [9j, however, all the 
solutions of had too large a mass to be the endpoint 
of the evolution of the unstable black string as conjec- 
tured in |icj | . Note however the recent work by Sorkin [I| 
(and related discussions in where it is discussed how 
for large enough spherical dimensions the static solutions 
found do indeed have lower mass, and could in principle 
be those conjectured by 0- 

In an attempt to resolve the issue of the fate of the 
linearly unstable black string, in Q a numerical dis- 
cretization of the Einstein equations in five-dimensions 
was presented. An SO (3) symmetry was assumed, effec- 
tively making the problem 2+1 dimensional, and hence 
one that could be solved on contemporary computer sys- 
tems. This code was employed to study the evolution 
of slight perturbations of a static, uniform black string. 
The simulation tracked the behavior of the system well 
into the non-linear regime, though unfortunately it could 
not be extended far enough to elucidate the final fate of 
the string. Nevertheless, the numerical solution revealed 
significant and rich dynamics, illustrating, in particular, 
that the string evolves toward a shape resembling large 
spherical black holes connected by thin strings. In terms 
of areal radius, the ratio of the maximum to minimum 
radius was about 10 by the end of the simulation. 

In this paper, we re-examine the data of Q, paying 
closer attention to the analysis of || and examine rele- 
vant features of the solution. In particular, we numer- 
ically find the expansion, shear and affine parameter of 
the black string horizon and discuss their behavior in con- 
nection to the conjecture of |8| where a scenario is ruled 
out as improbable, though as we argue here, it might very 
well be possible. The methods used are presented in sec- 
tion II, results in section III and conclusions in section 
IV. 
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II. METHODS 



Our starting point is to consider a spacetime described 
by the metric provided by the numerical simulation of . 
For simplicity we adopt the same gauge conditions used 
there (to avoid performing a numerical coordinate trans- 
formation, which in general would involve interpolating 
functions in space and time). The metric used in |lj takes 
the form 

ds 2 = {-a 2 +j AB (3 A f3 B )dt 2 + 2 lAB (3 A dx B dt 

+jABdx A dx B + jndil 2 (1) 

where x A = (r, z) and dD, 2 is the unit two-sphere met- 
ric. Our goal is to extract intrinsic properties of the 
event horizon and its generators. The horizon can be 
described by a surface r = R(t,z). To extract the 
sought-after information, the first step is to find the 
geodesic generators of the horizon. To this end, we con- 
struct the vector k a normal to the horizon defined as 
k a = {h,k r ,kz) = (— R, 1) — R')- Here an overdot de- 
notes differentiation with respect to t and a prime de- 
notes differentiation with respect to z. Since k a is null, 
it is tangent to the generators of the horizon; these are 
given by 



dz 
~dt 



r 



f3 r + R 



(2) 



Note however that these geodesies are not affinely param- 
eterized. To obtain the affinely parameterized geodesies 
one can exploit the fact that their tangent vectors l a obey 
the simple relation l a = e~ y k a where v is defined by 



dt f3 r + R 



(3) 



and h is given by 
1 d 



u 1 ( 
h= 2d~rK 



(4) 

Thus, the affine parameter A of these geodesies is found 
by integrating 



(5) 



dt fjr + R 



Next we turn our attention to obtaining the expan- 
sion and shear of these affinely parametrized generators. 
Since the spacetime is spherically symmetric, these quan- 
tities are given by just two components of V a h- Let x a 
be a unit vector in the two-sphere direction and let y a 
be a unit vector orthogonal to k a and to the two-sphere 
directions. Define the quantities A and B by 



A = 2x a x b VJ b 
B = -y a y b VJb 



(6) 
(7) 



which measure the distortion along spherical and longi- 
tudinal directions respectively. Then, the expansion and 
squared shear are given in terms of A and B by 



cr <?ab 



A-B 

\{A + 2B) 2 . 
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(8) 
(9) 



The equations above determine the affinely parametrized 
generators of the horizon and its expansion and shear. 
As we will see later, due to the exponential behavior of 
these quantities, it will turn out to be convenient to con- 
sider the logarithm of the affine parameter s = In A and 
rescaled quantities defined as follows 



{8,a ab ,A,B) = \(9,<j ab ,A,B) 
It then follows that 



ds 
It 



A-B 
\{A + 2Bf 



(3 r + R 
The quantity A is found by 

A = e { - s -^k a V a {\n ln ) 



(10) 

(11) 
(12) 

(13) 
(14) 



To find B we note that there is freedom to add a multiple 
of k a to y a and use that freedom to make y l vanish. It 
then follows that the components of y a are 



IT 



lz Z + 1 lrz R' + lrr {R'Y 



-1/2 



and y r = y z R 1 , from which 

B = ({y z ) 2 R" + k b T b ac yY 



(15) 



(16) 



Finally, we point out the following observation which 
will be crucial in the numerical evaluations to be carried 
out later. As we have seen above, 9 and <r ab (and their 
scaled counterparts) can be calculated independently. 
However they are related by Raychaudhuri's equation, 
which for the five-dimensional vacuum null case is 



dO 
dX 



¥ 2 



o, 



(17) 



while the corresponding equation for the rescaled quan- 
tities is 



dO 
ds 



(18) 



III. RESULTS 



The simulations performed in 1] provide a (partial) 
description of the spacetime describing a slightly per- 
turbed black string that evolves to a state where the 
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black string is quite deformed as judged by the appar- 
ent horizon. The surface that we use here as the event 
horizon is the boundary of the causal past of the region 
exterior to the apparent horizon at the latest time in the 
simulation (calculated using a generalization of the tech- 
nique described in yj), and it turns out that this surface 
is very well approximated by the apparent horizon for the 
entire length of the simulation. We initialize the affinc 
parameter to 1 (hence s = 0) at the beginning of the 
simulation. Among all generators, the one of particular 
interest is the one corresponding to the minimum radius, 
since it is there where any pinch off would occur. Due 
to the reflection symmetry of the initial data, the set of 
points on the horizon that at each time correspond to 
the minimal radius form a geodesic generator, and that 
generator stays at constant z. 

Figure^shows plots of s as a function of the simulation 
time coordinate t for the generator of minimum radius, 
from simulations with three different resolutions (to give 
some measure of the error in the results). Note that 
s has reached 50 near the end of the simulation. This 
corresponds to an affine parameter of about 10 21 . Thus, 
though the simulation has only taken a moderate amount 
of simulation time, it has gone for an enormous amount 
of affine parameter. 

Figure |21 shows plots of A and B vs s for the generator 
of minimum radius. Note that due to the rescaling this 
means that by the end of the simulation, the unrescaled 
expansion 9 and shear a a b are of order 10 -22 . Note too, 
that though the behavior of A and B as a function of s is 
quite dynamic throughout the simulation, since d/dX = 
X^d/ds it follows that viewed as functions of A even the 
rescaled expansion and shear would be regarded as slowly 
varying. 

We now turn to the behavior of the rescaled quanti- 
ties 9 and a a b- We can calculate the shear by simply 
evaluating (|12|) along the generator of minimum radius: 
Figure |3 shows the result of this as a function of s for the 
three simulations. In principle we can also calculate 9 in 
a similar manner using which says 9 is the differ- 

ence between A and B. However, the problem with this 
is that A and B are very close to one another in magni- 
tude, as can be seen from Figure |21 and so the difference 
will be dominated by numerical error. In fact, from the 
plot of A in Figure |21 and the assumed second order con- 
vergence of the numerical scheme, one can estimate that 
the error in A (and similarly for B), is on the order of 
5% at intermediate times (near the minimum of A) and 
grows to as much as 40% at late times. Therefore at late 
times the difference A — B will be completely swamped 
by numerical error. 

In order to alleviate this problem, we use Raychaud- 
huri's equation to integrate 9. However, in the form of 
(|18fl there is in general no stable direction of integration 
due to the linear and quadratic terms in 9; in other words, 
numerical errors in 9 could cause a blow-up integrating 
in either direction in s. Therefore, we eliminate the 9 2 
term using Qll[). and a ab a a b using iJHJ, though we keep 
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FIG. 1: s vs i for the horizon generator of minimum ra- 
dius, from three simulations of identical initial data though 
differing resolution (the curves are labelled by N r xN z , where 
N r (Nz) is the number of grid points in the r(z) direction). We 
show results from three different resolutions to demonstrate 
that we are in the convergent regime. Recall that the affine 
parameter A = e s , and therefore A is growing very rapidly 
with simulation time t. 

the linear term in 9 to give 

^.-0 + Ia 2 + B 2 -0. (19) 
ds 2 

In this form there is a stable direction of integration 
from large to small s, and the only 'sources' in A and B 
are quadratic and of the same sign, hence we are not as 
susceptible to numerical errors arising from the difference 
of similar numbers. One potential difficulty is in specify- 
ing initial conditions for the integration, which we need 
to do at the end of the simulation where 9 obtained from 
(lllll has the largest error. Fortunately the integration 
seems to be fairly insensitive to the initial condition in 9, 
as illustrated in Figure where we show the results of 
integration of (|19|l starting from both 9 calculated using 
(|lip. and 9 = (which is within the numerical margin of 
error at late times). 

IV. CONCLUSIONS 

We have found the behavior of the horizon of a black 
string in terms of the affine parameter, expansion and 
shear of the generators. The most significant aspects of 
our results are the extremely large values of the affine pa- 
rameter, and correspondingly the extremely small values 
of the (unrescaled) expansion and shear. In hindsight, 
this result is not too surprising as it is well known that 
for a Schwarzschild black hole there is an exponential re- 
lation between the Killing time coordinate and the affine 
parameter of the horizon generators. Since the straight 
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FIG. 2: The plot at the top depicts j4 as a function of s, while 
the plot at the bottom shows A and B vs. s from the highest 
resolution simulation, both for the horizon generator of mini- 
mum radius. What the latter plot demonstrates is that A and 
B are very similar in magnitude; in particular the difference 
A — B (which is just 9) at late times is completely dominated 
by numerical error, the magnitude of which can be estimated 
by using the data from the three simulations. Specifically, 
the error in A from the highest resolution simulation is quite 
small initially, is near 5% near the minimum of A, and has 
grown to around 40% by the end of the simulation (B has 
similar error). 



black string is just the Schwarzschild spacetime crossed 
with a circle, and since the spacetime treated here is ini- 
tially just a small perturbation of the straight string, it 
appears natural that very large values of the affine pa- 
rameter arc reached. The immediate conclusion is that 
in some sense the logarithm of the affine parameter is a 
more natural "dynamical time" than the affine parame- 
ter. 

What do our results say about the final fate of the un- 
stable black string? They certainly do not yet resolve the 
issue. The same three possible alternatives are still there: 
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FIG. 3: a ab a a t vs. s for the horizon generator of minimum 
radius. 



the string can (i) evolve to a static wiggly string, (ii) 
pinch off in finite affine parameter or (iii) pinch off in in- 
finite affine parameter. It is clear both from pj and from 
our work that the simulation has not gone far enough 
to read off the asymptotic behavior of the string. For 
that, a more robust simulation is needed. However, our 
work does shed extra light on the status of the alterna- 
tives, in particular that of alternative (iii). Recall that in 
Q it was correctly pointed out that alternative (iii) re- 
quires very small (unrescaled) expansion and shear and 
very slow dependence on affine parameter. Our results 
here show that precisely this is the sort of behavior that 
occurs, at least in some regime of the evolution of the 
black string. Moreover, this behavior is still in a regime 
of non-trivial dynamics. 

In fact, one might argue that alternative (iii) is the 
most plausible for lower spherical dimensions of the 
spacetime if the only non-uniform static solutions are 
those found in Q/. This fact would make alternative (i) 
unlikely as these solutions are too massive for spacetimes 
with dimension lower than 13. Additionally, option (ii) 
requires a singularity forming outside the horizon, and no 
hints of such behavior developing were seen in the simu- 
lations of [jj . Note that under alternative (iii) one might 
still have different scenarios. For instance, the 'minimal' 
circle might continue to shrink slowly with the trend ob- 
served at the late stages of the simulations of 1]. This 
requires the neighborhood of the string not to be de- 
scribable as a linearly perturbed uniform string since the 
ratio of length/area is above the critical length. On the 
other hand, if this is indeed the case, then the string may 
pinch off through a cascading sequence of instabilities on 
ever smaller scales, which again could happen in finite 
asymptotic time though infinite affine time. 

In a more quantitative vein, we can consider the possi- 
ble asymptotic behavior of the (rescaled) expansion and 
shear. Let p and c be positive constants. Then it is 
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lution. If p < 1 then this behavior results in pinch off 
in infinite affine parameter. If in addition, p > 1/2 then 
this pinch off occurs with finite J 9dX as assumed in 

If such a power law is the asymptotic behavior of the 
unstable black string, then the simulations have not yet 
run for long enough to be sufficiently far into the asymp- 
totic regime to try to calculate the exponents. Nonethe- 
less, in Figures [3] and 0] the quantities a ab a a b and 9 peak 
and then decrease. This decrease may be the beginning 
of a power law tail to be found in an improved black 
string simulation. 
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FIG. 4: 6 vs. s for the horizon generator of minimum radius, 
calculated by integrating (1191 . For the integrations shown on 
the top figure, the initial value of (which is at the largest 
s, as we integrate from large to small s) is 0, while on the 
bottom figure the initial 9 was calculated using [11 1 1 . We show 
curves calculated using two different sets of initial conditions 
to demonstrate that the integration is fairly insensitive to the 
initial value of 6; hence, expect near the end of the simulation, 
the error in the integrated 9 is comparable to that of 5- ab a a b 
shown in Figure [3] — specifically around 10% near the maxima 
for the highest resolution simulation, assuming second order 
convergence. 



consistent with equation 1|18|) for both 9 and a ah <j a \, to 
asymptotically approach c 2 s~ 2p . This would result from 
both A and B approaching —cs~ p asymptotically, but 
their difference approaching c 2 s~ 2p . That this sort of 
behavior is possible was recognized in Q but there it 
was thought that such behavior was somehow unnatu- 
rally slow and indicative of an approximately static so- 
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